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ABSTRACT 

We present new spatially resolved observations of the dust thermal emission at 7 mm from the Fomalhaut debris disk obtained with 
the Australia Telescope Compact Array. These observations provide the longest wavelength detection of the Fomalhaut debris disk 
to date. We combined the new data to literature sub-mm data to investigate the spectral index of the dust thermal emission in the 
sub-millimeter and constrained the g-slope of the power-law grain size distribution. We derived a value for q = 3.48 ± 0.14 for grains 
with sizes around 1 mm. This is consistent with the classical prediction for a collisional cascade at the steady-state. The same value 
cannot be explained by more recent collisional models of planetesimals in which either the velocity distribution of the large bodies or 
their tensile strength is a strong function of the body size. 
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1. Introduction 

Dusty disks around nearby main-sequence stars have been 
known for almost three decades (Aumann et al. 1984). These are 
typically called debris disks, because the detected dust is thought 
to be produced by collisions between km-sized planetesimals, 
leftovers of the planetary formation process (see Wyatt 2008). 

The highly destructive mutual encounters between planetesi- 
mals, triggered either by dynamical interaction with one or more 
planets in the system (Mustill & Wyatt 2009) or by self-stirring 
(Wyatt 2008, and references therein), produce objects down to 
very small sub-/im dust grains. The size distribution of the solids 
produced by these collisions provides information on some of 
the main properties of the invisible planetesimals. For example, 
collisional models of debris disks provide different predictions 
for the grain size distribution depending on the dynamical state 
(velocity dispersion, e.g. Pan & Schlichting 2011) or physical 
conditions (e.g. tensile strength, Durda & Dermott 1997) of the 
large bodies. 

Whereas the jjm sized (and smaller) grains strongly inter- 
act with the stellar radiation via radiation pressure, the dynam- 
ics of ~ 0. 1 millimeter sized (and larger) pebbles are dominated 
by collisions and gravitational interactions. For this reason these 
pebbles remain located with the planetesimals (e.g. Wyatt 2005, 
Corder et al. 2009, Wilner et al. 201 1). The spectral energy dis- 
tribution (SED) of debris disks in the millimeter range can be 
used to constrain the size distribution of these pebbles and pro- 
vide a powerful tool to test the collisional models of planetesi- 
mals (e.g. Vandenbussche et al. 2010). 

In this paper, we present new spatially resolved 7 mm 
observations of the Fomalhaut debris disk obtained with the 
Australia Telescope Compact Array (ATCA). The Fomalhaut 



main-sequence star (a PsA, spectral type A4) is surrounded 
by one of the closest known dusty debris disks (distance of 
7.66±0.04 p c from Hipparco s, van Leeuwen 2007 1. Its age 
(~ 200 Myr, |Di Folco et al.||2004| i greatly exceeds the typical 
timescale for the survival of pristine circumstellar dust grains (< 
5-10 Myr) because of interaction with gas in primordial disks 
(see Pinilla et al. 2012, and Williams & Cieza 201 1 for a recent 
review), so that a continuous replenishment of dust through col- 
lisions of planetesimals is needed. 

The geometry of the Fomalhaut debris disk is characterized 
by a very sharp inner edge at a radial distance of about 133 AU 
from the central star; the geometric center of the disk is also off- 
set by about 15 AU from Fomalhaut (Kalas 2005); both these 
features were soon recognized as likely signatures for one or 
more planets shaping the morphology of the disk. A few years 
later, a planetary mass object, Fomalhaut b, has been imaged by 
Kalas et al. (2008) using the Hubble Space Telescope Advanced 
Camera for Surveys coronograph. This disk is therefore a bench- 
mark for investigating the interaction between a planetary sys- 
tem and its debris disk (Chiang et al. 2009). 

The dust continuum emission of this source has been mea- 
sured in the sub-mm at about 0.35 mm with the CSO/SHARC 
II telescope (Marsh et al. 2005), at 0.45 mm and 0.85 mm with 
JCMT/SCUBA (Holland et al. 1998, 2003), and at 1.3 mm with 
SEST (Chini et al. 1991). The fluxes measured in the sub-mm 
make the Fomalhaut system one of the brightest debris disks at 
these long wavelengths. Therefore this is one of the very few 
cases where current technology allows an investigation of the 
SED to wavelengths longer than ~ 1 mm. 

Together with the ATCA observations of the ft Pic debris 
disk published in Wilson et al. (2011), our new observations 
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provide the longest wavelength detection of any debris disk to 
date. This gives us the advantage of probing grains as large as 
about 1 cm, which are virtually unaffected by the stellar radia- 
tion. At the same time, the addition of a flux measurement at the 
long-wavelength end of the SED allows us to better probe the 
Rayleigh-Jeans tail of the spectrum, thus minimizing the effect 
of any possible bias introduced by modeling the dust tempera- 
tures from the observed data. The new 7 mm data also provide a 
better lever arm in wavelength that mitigates the impact of abso- 
lute calibration uncertainties on the measurement of the spectral 
index. We present an analysis of the long-wave emission of the 
Fomalhaut debris disk and discuss its impact for the theories of 
collisional systems of planetesimals. 

2. Observations 

We observed the Fomalhaut debris disk at 7 mm with ATCA 
and its Compact Array Broadband Backend (CABB) digital filter 
bank (Wilson et al. 2011). The observations were carried out 
on September 15, 16, and 17, 2011. The ATCA array was in 
the H75 configuration, which provides baselines between 3 1 and 
89 m. To obtain the best sensitivity in the continuum we set the 
correlator to cover the full 4-GHz CABB effective bandwidth. 
The specific mean wavelength of the observations is 6.66 mm. 

The raw visibilities for each day were calibrated and edited 
with the MIRIAD software package. The bandpass and absolute 
flux calibrations were made with 1921-293 and Uranus, respec- 
tively. The complex gain calibration was performed on the cal- 
ibrator 2255-282, and the Fomalhaut debris disk was typically 
observed for 10 minutes followed by a 2 minute integration on 
the gain calibrator. The total integration time on the Fomalhaut 
debris disk was about 12 hours. On the three days of observa- 
tions, the measured fluxes of the 2255-282 calibrator were found 
to be consistent within about 4% from the mean value (3.61 Jy). 
To account for the systematic effects that typically affect interfer- 
ometric observations (e.g. baseline errors, pointing errors, etc.), 
we considered an uncertainty of about 15% on the calibrated flux 
at ~ 7 mm. The l<x noise level in the final map is 23 yuJy. The 
angular resolution of our image is about 14.3 x 10.7" (PA m - 
76 deg), and the primary beam full- width-at-half -power (FWHP) 
is about about 65" centered at the location of the Fomalhaut star 
(Figure [T}. 

3. Results 

3.1. ATCA map at 7mm 

Figure [T] shows the new map of the Fomalhaut disk obtained 
with ATCA at 7 mm. The disk is clearly resolved into two main 
components, or lobes, northwest (NW) and southeast (SE) from 
the central star, respectively. These two lobes are detected at a 
signal-to-noise ratio of ~ 6. The emmision from the SE lobe is 
unresolved at the angular resolution of the ATCA observations. 

The two-lobe morphology seen in the ATCA map is 
very similar to the one observed at 0.85 mm by Holland et 
al. (1998,2003) with the JCMT/SCUBA single-dish telescope at 
about the same angular resolution (« 14"). The position angle 
of the two lobes, evaluated at the peak of the surface brightness, 
is also nearly identical (« 162, 159 degrees from the ATCA and 
SCUBA maps, respectively). The angular resolution of our ob- 
servations is not high enough to resolve the ringlike structure as 
observed by Marsh et al. (2005) at 0.35 mm with the SHARC II 
camera at the CSO single-dish telescope. 

In contrast to imaging from single-dish telescopes, interfer- 
ometric imaging can filter out a fraction of the source emission. 
The largest angular scales that can be probed by interferometric 



observations are determined by the shortest spacings between 
any pair of antennas in the array. The ATCA observations pre- 
sented here are sensitive to emission at angular scales up to ~ 
50". Since this is very close to the largest scales for the emission 
detected for the Fomalhaut disk by single-dish telescopes in the 
sub-mm, no significant filtering of emission from the disk is ex- 
pected in our ATCA image. We verified this by taking a simple 
narrow ring model similar to that used by Marsh et al. (2005) for 
the Fomalhaut disk; we simulated interferometric observations 
with the same coverage of the (u, v)-plane as for our ATCA ob- 
servations, and we found that more than 99% of the flux from 
the model was recovered by the simulated observations. 

In Figure [T]the blue cross and orange star symbols show the 
locations of the disk center (the halfway point between the peaks 
of the NW and SE lobes) and of the star, respectively; the dis- 
crepancy between the positions of the disk center and of the star, 
i.e. about 2", or 15 AU, is consistent with that found by Kalas 
et al. (2005) with HST The fact that this same discrepancy is 
also seen at millimeter wavelengths confirms that the asymme- 
try involves the planetesimal ring, rather than just the smallest 
fim sized grains seen in the optical. 

The NW lobe appears to be more extended, showing two 
possible asymmetric structures toward east and south. These 
possible asymmetries are observed at low signal-to-noise ratios 
(« 2 - 4), and more sensitive observations are needed to deter- 
mine the real nature of these apparent structures. The 4cr-contour 
in the southern asymmetry of the NW lobe extends to about 4" 
from the location of the Fomalhaut star (see Figure[T]i; this emis- 
sion, if real, might be caused by an ionized gas component (e.g. 
from a stellar wind or stellar corona) rather than by dust thermal 
emission from the debris disk; sensitive observations at centime- 
ter wavelengths are needed to measure the spectral index of this 
emission and determine its nature. 

The integrated flux, including the central emission, is about 
400 ± 64 fiiy (the uncertainty accounts for the absolute flux un- 
certainty), whereas the central emission is 92 //Jy, or 23% of 
the total emission. The possible effect of the uncertainty on the 
naure of the central emission on our analysis will be discussed 
in Section l3~2l 

3.2. Millimeter spectral index and slope of the grain size 
distribution 

The angular resolution of the new ATCA observations does 
not allow a detailed investigation of the morphology of the 
Fomalhaut debris disk. At the same time, the combination of the 
integrated flux derived at 7 mm with the sub-mm literature data 
allows us to probe the Rayleigh-Jeans tail of the dust thermal 
emission and derive information on the size distribution of the 
emitting grains. 

Thermal dust emission from debris disks is optically thin, 
and therefore the flux density F v oc fi v (7d ust ) x k v x M^ ust /d 2 , 
where fi v (7d U st) is the Planck function at the characteristic dust 
temperature Tdust, « v is the dust opacity coefficient, Md us t is 
the dust mass, and d is the distance. At long (sub-)millimeter 
wavelengths the dust opacity k v oc yP, and the Planck function 
B v (7 'dust) x v a " with a P i x 2. Hence, a mm ~ ffpi+yS, where 
the (sub-)mm spectral index of the millimeter SED (F v oc y amm ). 

The yS-parameter is related to the slope q of the power-law 
differential size distribution of the dust grains (dn(a) oc cC q da). 
Draine (2006) derived the analytical formula /? « (q — 3) • f3 s , 
where /3 S is the dust opacity spectral index in the small par- 
ticle limit. This is valid for 3 < q < 4, and as long as the 
size distribution of grains is a single power-law over a broad 
enough interval of grain sizes around the observed wavelengths, 
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Fig. 1. ATCA image at ~ 7 mm of the Fomalhaut debris disk corrected for the attenuation of the primary beam pattern up to the half 
power point (« 33" from the pointing center). White contours are drawn at (-2, 2, 4, 6) x <x, where lcr ss 23 yuJy is the rms noise 
level on the map. The orange star symbol identifies the location of the Fomalhaut star (a, 5)=(22:57:39. 05, -29:37:20.05), whereas 
the disk center is indicated by the blue cross, whose size corresponds to the angular resolution of the AT CA observations (=a 2"). 
Coordinates in the map have been corrected for the Fomalhaut proper motion (329/-164 mas yr _1 , from van Leeuwen |2007| , and 



are thus correct for the 2000.0 epoch. The red ellipse in the bottom left corner denotes the synthesized beam of the observations 
(FWHM * 14.3x10.7", PA * -76 deg). 



i.e. a m ; n << 1 mm, a max >> 1 mm (Draine 2006). This relation 
has been found to be very accurate for different models of dust 
considered in the literature to interpret the mm-wave emission of 
young circumstellar disks (see e.g. D'Alessio et al. 2001, Ricci 
et al. 2010a,b,201 1). Combining the relations written above, one 
can obtain a simple relation linking the slope q of the grain size 
distribution to the spectral indices a^, an, and f} s \ 



■ an 



+ 3, 



(1) 



which is valid for 3 < q < 4. 

Assuming that particles in debris disks are composed of sim- 
ilar materials as interstellar grains and as are expected in primor- 
dial disks, the dust opacity spectral index in the small particle 
limit p s = 1.8 ± 0.2, where the uncertainty accounts for the dif- 
ferent /^-values measured in different diffuse interstellar clouds 
and in dense molecular clouds (see Draine 2006 and references 
therein). The spectral index of the Planck function an depends 
on the temperature of the emitting dust, which is equal to 2 in 
the Rayleigh-Jeans limit. Holland et al. (2003) derived an in- 
terval of w 40 - 50 K for the temperature of the large (>1 mm) 
grains observed in the millimeter, by fitting infrared and sub-mm 
data. This is consistent with the best-fit value of 42 K derived by 
Marsh et al. (2005). At these temperatures, the spectral index of 
the Planck function, calculated between 0.35 and 6.66 mm, is 
an = 1-84 ± 0.02. Once the f} s and an parameters are estimated, 
Equation[T]allows us to constrain the slope q of the grain size dis- 
tribution by measuring the spectral index ff„ of the (sub-)mm 
SED. 

Figure [2] shows the SED of the Fomalhaut debris disk be- 
tween about 0.35 mm and 6.66 mm. Here we consider the flux 
measured at 1.3 mm by Chini et al. (1991) only as a lower 
limit, because their SEST observations have been performed 
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Fig. 2. Sub-millimeter SED of the Fomalhaut debris disk. Data 
at 0.35, 0.45, 0.85 are from Marsh et al. (2005) and Holland et 
al. (2003); the plotted uncertainties for these data are 20%, 20%, 
and 15% of the measured flux, respectively; the data point at 
1.3 mm from Chini et al. (1991) is taken as a lower limit and is 



not considered for the fit (Section 3.2 1; the blue data-point indi- 
cates the integrated flux density at 6.66 mm from this paper. The 
value of the spectral index measured between 0.35 and 6.66 mm 
is also shown in blue. 



with a single beam (FWHP =s 24") pointing at the stellar po- 
sition. Consequently, a significant fraction of the disk emission 
has been missed. Gaspar et al. (201 1) instead included the Chini 
et al. (1991) detection in their analysis which explains the rel- 
atively high value of the spectral index they obtained, between 
0.35 and 1.3 mm (a: 3.1 + 0.3). In Figure |2]we also report the 
measured value of the spectral index ao.35-6.66mm = 2.70 + 0.17. 
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Inserting this value in Equation [T] yields a value for the slope of 
the grain size distribution q = 3.48 ± 0.14, where the uncertainty 
is derived by propagating the uncertainties on the a mm ,ap\,j3 s 
parameters through Equation[T] The central emission at ~ 7 mm 
discussed above is associated with ionized gas and not with the 
debris disk, the constrained value of q - 3.53 + 0.15 would be 
fully consistent with the value reported above (the correspond- 
ing spectral index ao.35-6.66mm = 2.80 ± 0.20). The effect of the 
central emission on the derived value of q is therefore negligible. 
We also checked that our result is not significantly biased by our 
choice of the interval in wavelength, e.g. due to deviations from 
the Rayleigh-Jeans regime especially at ~ 0.35 - 0.45 mm. By 
repeating the same analysis in the 0.85-6.66 mm interval, we ob- 
tained q = 3.43 + 0.17, which is consistent with the value derived 
between 0.35 and 6.66 mm, although with a larger error because 
of the shorter lever arm in wavelength. 

4. Discussion 

Different collisional models of debris disks have provided pre- 
dictions for the grain size distribution. In general, the destruc- 
tive mutual encounters between planetesimals produce objects 
down to very small sub-yum grains; also, for most of the mod- 
els, the predicted size distribution is satisfactorily approximated 
by a power-law with a slope q assuming values between 3 and 
4. Therefore, the assumptions used to derive Equation [T] are sat- 
isfied and the constraints we derived in the last section for the 
g-slope of the Fomalhaut disk can be used to test the predictions 
of different models. 

This method allows one to determine the slope of the grain 
size distribution only in a relatively narrow interval of sizes be- 
tween ~ 0.1 mm and ~ 1 cm, since these grains are the most 
efficient emitters in the sub-millimeter. Solids much larger than 
~ 1 cm are not probed by our observations. But our method has 
the advantage of tracing grains that are virtually unaffected by 
the stellar radiation. Therefore, this provides more robust con- 
straints for the grain size distribution than is currently possible 
by analyzing the infrared emission, which probes the smaller and 
more pressure radiation affected /vm-sized particles. 

The value of the g-index derived in the last section, i.e. 
q = 3.48 ± 0.14, is consistent with the classical prediction of 
q = 3.51 for a collisional cascade at the steady-state, as a re- 
sult of the collective dynamical interaction of particles caused 
by inelastic collisions and fragmentation (Dohnanyi 1969). The 
standard theoretical treatment of collisional cascades assumes a 
single constant tensile strength^ and velocity dispersion for all 
bodies regardless of size. Both these quantities, and especially 
their possible variation with the size of the solids, are known to 
play a role in the predicted slope of the grain size distribution. 

Recently, Gaspar et al. (2011) used numerical models of 
collision-dominated debris disks to investigate the effects of 
model parameters on the evolution of the grain size distribution. 
These authors showed that the slope of the tensile strength curve 
is the parameter that affects the slope of the grain size distribu- 
tion the most (see also Durda & Dermott 1997). In particular, 
the authors found that at the age of the Fomalhaut disk, values 
of q > 3.82 are obtained if the strength curve varies with the 
solids size as aT x or steeper. Adopting Gaspar et al.'s fiducial 
model for the Fomalhaut disk, our observations rule out such 
steep tensile strength curves at the >95% confidence level. This 
discrepancy can be partially alleviated if other degenerate model 
parameters are simultaneously varied, but our observations al- 

1 The tensile strength of a body is defined as the minumum energy 
per unit mass required to disrupt it catastrophycally. 



ways tend to favor shallower strength curves. In particular, the 
slope q obtained by our analysis is consistent with the prediction 
of q — 3.65 of the Gaspar et al. (2011) fiducial model, where 
they adopted a tensile strength curve varying with body size as 
a" 038 . 

Pan & Schlichting (2011) have recently extended the classi- 
cal treatment of collisional cascades (Donhanyi 1969) by relax- 
ing the hypothesis of a single constant velocity dispersion, and 
solving self-consistently for the body size and size-dependent 
velocity distributions at the steady-state. They derived a veloc- 
ity distribution v(a) cc a 5 and a correspondingly very steep 
slope of the body's size-distribution q = 4 in the sizes of inter- 
est for the comparison with our results for the Fomalhaut disk. 
This predicted value of q is ruled out by our analysis at the 99% 
confidence level. Our analysis therefore suggests that the vari- 
ation in body size of the velocity distribution of planetesimals 
in the Fomalhaut debris disk is weaker than predicted by the 
model of Pan & Schlichting (201 1). This might be explained by 
the gravitational stirring expected in the Fomalhaut planetesimal 
belt caused by the known massive exo-planet (Fomalhaut b) de- 
tected close to the debris disk (Mustill & Wyatt 2009). 

The application of this method to other debris disks with and 
without massive exo-planets will be extremely interesting for in- 
vestigating possible signs of interactions between planets and 
planetesimal belts. This will be made possible in the near future 
by the very high sensitivities and angular resolutions offered by 
ALMA and the EVLA at mm wavelengths. 
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